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THE  DERIVATION  OF  THE  MAPPING  EQUATIONS  AND  DISTORTION 


FORMULAE  FOR  THE  SATELLITE  TRACKING  MAP  PROJECTIONS 

by 

Gregory  Cote  Arnold 
(ABSTRACT) 

Rigorous  derivations  of  the  Satellite  Tracking  conic  and  cylindrical 
projections  are  presented.  The  first  fundamental  quantities  of  the 
projection  surface  parameters  as  functions  of  spherical  earth  para¬ 
meters  are  developed.  Using  these  quantities,  newly  derived  in  this 

paper,  general  equations  for  the  distortion  in  length,  area,  and 
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azimuth  are  developed.  Examples  of  the  graticule  and  distortion 
values  are  given  for  the  Landsat  orbit. 
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1 .  INTRODUCTION 


The  growing  use  of  remotely  sensed  earth  imagery  has  prompted 
interest  in  a  new  type  of  map  projection,  one  in  which  some  relation¬ 
ship  between  satellite  position  and  the  earth  is  preserved  during  the 
transformation,  rather  than  the  preservation  of  some  feature  or  prop¬ 
erty  of  the  earth  itself.  A  good  example  is  the  Space  Oblique 
Mercator  described  by  Snyder  (1978)  on  which  conformality  is  preserved 
along  a  satellite  groundtrack. 

The  Satellite  Tracking  projections  have  the  property  of  portraying 
all  groundtracks  of  a  particular  satellite  as  straight  lines.  The 
advantage  of  this  is  that,  once  the  graticule  is  constructed,  knowing 
the  location  on  the  sphere  of  any  subpoint  allows  the  entire  revolu¬ 
tion  to  be  drawn. 

1 . 1  Background 

In  his  article  "Map  Projections  for  Satellite  Tracking",  Snyder  (1981) 
introduces  the  cylindrical  and  conic  Satellite  Tracking  projections. 
The  article  briefly  discusses  the  concept  of  satellite  apparent 
longitude  from  which  the  relationships  between  the  revolving  satellite 
and  the  rotating  earth  are  derived,  and  includes  the  development  of 
the  mapping  equations  and  examples  of  the  graticules  of  both  projec¬ 
tions. 
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Although  the  mapping  equations  are  developed,  the  complete  relation¬ 
ship  between  coordinates  on  the  sphere  and  coordinates  on  the  projection 
surface  is  not.  To  define  this  relationship,  and  subsequently  derive 
general  formulae  for  expressing  various  types  of  distortion  on  the 
projection  surface,  the  fundamental  quantities  of  each  projection  must 
be  known.  These  are  not  determined  by  Snyder. 

1.2  Objective 

The  objective  of  this  thesis  is  the  derivation  of  general  distortion 
formulae  for  the  Satellite  Tracking  map  projections.  This  is  accom¬ 
plished  by  first  deriving  the  fundamental  quantities  of  the  projections 
which  relate  the  change  in  on  the  sphere  t«  the  change  in  the  pro¬ 
jection  parameters  on  the  plane.  From  quantities,  general  equa¬ 

tions  for  length,  area,  and  angular  distortion  are  derived.  The  equa¬ 
tions  used  in  the  computation  of  the  fundamental  quantities  and  distor¬ 
tion  formulae  are  given  in  Pearson  (1977). 

1.3  Scope 

After  a  brief  explanation  of  terms  and  concepts  as  given  in  Richardus 
and  Adler  (1972),  a  complete  rederivation  of  the  two  projections  is 
presented.  The  derivations  are  included  here  for  two  reasons:  first, 
each  step  is  completely  explained,  which  should  make  them  easier  to 
follow  than  the  original;  and  second,  the  resulting  mapping  equations 
are  the  basis  for  the  determination  of  the  distortion  formulae. 

In  the  chapter  on  distortion  analysis,  the  fundamental  quantities  and 


distortion  equations  are  derived.  Finally,  numerical-  examples  using 
the  Landsat  orbit  are  given. 

Changes  in  the  definition  of  some  variables  will  affect  the  appearance 
of  several  equations.  Specifically,  motion  of  the  satellite  along  its 
orbit  will  be  calculated  from  the  ascending  node,  the  customary  point 
of  reference,  rather  than  the  descending  node.  Also,  Snyder  refers  to 
a  non-rotating  earth  which  he  later  allows  to  rotate  to  complete  his 
description  of  satellite  motion.  Here,  the  non-rotating  earth  concept 
is  replaced  by  referencing  the  satellite  to  an  astronomical  coordinate 
system. 

Three  conditions  are  assumed  in  the  derivation  of  the  mapping  equations 

1.  The  earth  is  spherical  and  of  uniform  mass, 

2.  the  orbit  is  circular,  i.e.,  satellite  velocity  is  constant, 

3.  the  rotation  axes  of  the  selected  astronomical  and  earth-fixed 
coordinate  systems  are  coincident. 

Plotting  the  groundtrack  is  equivalent  to  orbit  determination  where 
the  above  conditions  serve  as  the  math  model  of  the  orbit.  Recogni¬ 
zing  the  limitations  of  such  a  crude  model,  extrapolation  on  any 
groundtrack  should  be  restricted  to  the  revolution  containing  the 
known  satellite  subpoint. 

1.4  Satellite  Apparent  Longitude 

Two  types  of  longitude  are  used  in  the  derivation  of  the  mapping 
equations:  first,  earth-fixed,  or  geodetic,  longitude  (A),  and  second, 
satellite  apparent  longitude  (A"),  defined  in  section  3.1.3.  The 


meridians  on  the  projection  represent  geodetic  longitude.  However,  in 
order  to  compute  the  parallel  spacing  necessary  to  allow  straight  line 
groundtracks,  the  motion  of  the  satellite  must  be  taken  into  consider¬ 
ation.  Therefore,  the  derivation  of  the  mapping  equations  for  41  will 
be  based  on  A'  instead  of  A. 


2.  MAP  PROJECTION  CONCEPTS 


2.1  The  Projections 

The  transformations  are  applied  to  positions  on  a  spherical  earth  of 
radius  R,  defined  in  the  latitude,  longitude  ($,A)  coordinate  system. 

The  two  projection  surfaces,  the  cylinder  and  the  cone,  are  oriented 
relative  to  the  sphere  in  the  normal  aspect,  i.e.,  their  axes  coincide 
with  the  polar  axis  of  the  earth.  The  plane  coordinate  systems  are 
chosen  such  that  the  parametric  lines  of  that  system  correspond  with 
the  projection  of  the  parametric  lines  of  the  sphere.  On  the  cylinder, 
where  the  graticule  is  composed  of  perpendicular  straight  lines,  the 
Cartesian  (x,y)  system  is  used.  On  the  cone,  where  meridians  are 
straight  lines  converging  at  its  apex,  and  parallels  are  arcs  of  concen¬ 
tric  circles,  the  (p,0)  system  is  used.  The  line  of  contact  between 
the  projection  surface  and  the  sphere  is  a  parallel  of  latitude  along 
which  scale  is  true.  A  secant  cone  or  cylinder  will  have  two  standard 
parallels,  and  a  tangent  surface  will  have  one. 

2.2  Distortion  On  Map  Projections 

Regardless  of  the  type  of  transformation,  all  of  the  relationships 
existing  on  the  sphere  cannot  be  duplicated  on  the  plane.  Any  pro¬ 
jection  will  contain  distortions  in  distance,  direction,  size  or 
shape.  The  Satellite  Tracking  projections  sacrifice  some  desirable 
properties  in  order  that  satellite  groundtracks  are  preserved  as 


straight  lines. 
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2.2.1  Length  distortion 

Length,  or  scale,  distortion  is  a  measure  of  the  difference  in  ler 
of  a  line  on  the  projection  surface  compared  with  the  length  of  tl 
same  line  on  a  sphere  which  has  been  reduced  to  map  scale.  It  is 
usually  expressed  as  a  scale  factor,  the  ratio  of  the  two  distance 

Scale  factor  =  distance  on  the  plane 
distance  on  the  sphere 

The  scale  factor  along  a  standard  parallel  or  meridian,  where  seal 
true,  is  1.  A  scale  factor  greater  than  1  indicates  an  increase  i 
length  on  the  projection;  a  value  less  than  1  indicates  a  decrease 
in  length. 

On  most  projections,  the  scale  factor  is  dependent  on  the  azimuth 
line.  Formulae  to  be  developed  here  are  for  distortion  along  the 
meridians  (Mm)  and  along  the  parallels  (Mp). 

Conformality  is  the  property  which  maintains  the  shapes  of  differe 
tially  small  areas  on  a  projection.  This  means  that,  although  the 
magnitude  of  the  scale  distortion  changes  across  the  map,  at  any  p 
Mm  =  Mp. 

This  condition  will  be  used  in  the  derivation  of  the  mapping  equat 

2.2.2  Area  distortion 

Area  distortion  (Da),  also  expressed  as  a  ratio,  compares  an  area 

the  map  with  the  same  area  on  the  earth  at  map  scale. 

Da  =  area  on  the  plane 
area  on  the  sphere 
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Because  the  parametric  curves  of  the  two  plane  coordinate  systems  are 
orthogonal,  the  area  bounded  by  two  differentially  close  pairs  of 
curves  is  the  product  of  the  lengths  of  the  lines  separating  them. 

Da  =  Min  Mp 

2.2.3  Angular  distortion 

Angular  distortion  is  measured  by  comparing  an  azimuth  on  the  sphere  to 
its  projected  azimuth  on  the  plane.  For  the  Satellite  Tracking  pro¬ 
jections,  it  is  a  function  of  azimuth  and  latitude.  Angular  distortion 
will  be  zero,  i.e.,  sphere  azimuth  equals  plane  azimuth,  wherever  the 
condition  of  conformality  is  satisfied;  otherwise,  some  distortion  will 
be  present. 

At  each  point  on  the  sphere  there  exists  a  pair  of  perpendicular  lines 
that  will  remain  perpendicular  during  the  transformation;  along  these 
lines,  angular  distortion  is  zero. 


3.  THE  DEVELOPMENT  OF  THE  MAPPING  EQUATIONS 


3.1  Elements  Of  The  Satellite  Orbit 

Unlike  projections  which  deal  only  with  parameters  of  the  earth-fixed 
coordinate  system,  the  Satellite  Tracking  projections  must  take  into 
account  the  motions  of  satellite  revolution  and  earth  rotation.  This 
is  because,  under  the  above  assumptions,  the  satellite  orbital  plane  is 
stationary  with  respect  to  the  stars,  with  its  earth  orientation  con¬ 
stantly  changing  due  to  earth  rotation.  Because  the  transformation 
between  these  two  systems  is  a  rotation  about  the  Dolar  axis,  only 
longitude  is  affected. 

To  develop  the  projections,  satellite  motion  along  its  orbit  must  be 
defined  in  terms  of  earth-fixed  coordinates  (<p,X).  This  is  done  by 
first  defining  the  orbit  in  astronomical  coordinates  ($,A),  and  then 
applying  the  transformation  to  (4>,A).  The  elements  of  the  orbit  are 
shown  in  Figure  1. 

3.1.1  Distance  along  the  satellite  track  as  a  function  of  <j 

Applying  the  sine  law  to  triangle  ACD  in  Figure  1  gives  the  relation¬ 
ship  between  the  angular  distance  traveled  by  the  satellite  (a)  as  a 
function  of  $  and  the  inclination  of  the  orbital  plane  (t). 

sina  =  sin<t>  (1) 

sini 

Under  the  conditions  imposed  for  the  derivation,  astronomical  latitude 
equals  geodetic  latitude.  Therefore,  4>  =  $,  and  equation  (1)  mav  be 
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rewritten  as 


sina  =  sin  6  (2) 

sin  i 

3.1.2  Distance  along  the  satellite  track  as  a  function  of  A 

Applying  the  sine  law  to  triangle  ABC,  and  the  cosine  law  to  triangle 
ACD  in  Figure  1  provide  the  following  relationships.  In  both  cases,  A 
is  measured  from  the  ascending  node. 


sina  = 

cos  $ 

sinA 

(3) 

cos  \ 

cosa  ■ 

cos$ 

cosA 

(4) 

Dividing  equation 

(3)  by  (4)  eliminates  $  and  expresses  a  as  a 

function  of  A  only 

• 

tana  =» 

tanA 

(5) 

COS  1 

3.1.3  Satellite  apparent  longitude 

The  geodetic  longitude  of  a  point  beneath  the  satellite  orbit  is 
equal  to  the  astronomical  longitude  plus  the  combined  effects  of  sat¬ 
ellite  and  earth  motion  on  earth-fixed  longitude  (AX). 

X'  =»  A  +  AX  (6) 

where  the  primed  X  indicates  longitude  of  the  satellite,  known  as 
satellite  apparent  longitude.  The  angular  distance  along  the  orbit 
covered  over  time  t  is 

a  -  2irt  (7) 

P 


Jb. 
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where  p  is  the  satellite  period  of  revolution.  During  the  same 
interval,  the  earth-fixed  longitude  of  the  satellite  subpoint  changes  by 

AX  =  -2irt  (8) 

P 

where  P  is  the  earth  period  of  rotation.  AX  is  negative  because  the 
earth  rotates  in  the  positive  X  direction,  therefore,  the  longitude  of 
the  point  is  always  decreasing.  Equate  (7)  and  (8)  to  obtain 

AX  =  -  £  a  (9) 

P 

Substitute  equation  (9)  into  (6)  to  obtain  the  final  form  of 
satellite  apparent  longitude 

X'  =  A  -  £  a  (10) 

P 

3.1.4  The  relationship  between  X*  and  <j> 

Groundtracks  on  these  projections  will  be  constrained  as  straight  lines 
by  spacing  the  parallels  as  a  linear  function  of  the  satellite  apparent 
longitude  as  computed  in  equation  (10).  To  do  this,  dX'/di j>  is  required. 
Differentiating  equations  (10),  (2),  and  (5) 

dX'  =  dA  -  p  dq  (11) 

d<j>  d$  P  d$ 

da  *  cos4>  (12) 

d$  sini  cosa 

dA  =  sec2q  cos\  da  (13) 

d$  sec2A  d$ 


•  < 


•  t 


»  « 


•  « 


I  « 


»  ( 


I  < 


I  < 


l 


< 
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Substituting  (12)  and  (13)  into  (11) 


dX^  =  /sec2q  cos;  “  P  A  cosi$ 
d<J>  sec2A  P  J  sin;  cosa 

dX'  =  /  cosi _  -  j>\ _ c 

d$  lcos2a  (1  +  tan2A)  P J  sin; 


Reducing  terms  in  the  last  equation: 


cos2a  (1  +  tan2A)  =  cos2a  +  cos2a  tan2A 

=  cos2a  +  cos2a  tan2a  cos2i 
=  1  -  sin2a  sin2i 


1  -  sln2<fr  sin" 


from  (5) 


from  (2) 


sinzi 


COS2(j) 


sini  cosa  ■  sin 


sini/l  -  sin2 (A1* 

l  TirfJ 

(sin2i  -  cos2^)1* 
(cos2#  -  cos2 ;  )^ 


from  (2) 


Continuing  the  derivation: 


dX  '  =  / co si.  - 
d <t>  \cos2$  ? J  ( 


cos^A  -  COSM, 


dX'  ■  cos;  -  (p/P)  cos2<i 


d<()  cos$  (cos2ij)  -  cos2;) 


Given  dX'  and  d$  on  the  sphere,  the  azimuth  (A)  of  a  great  circle  at 
latitude  $  is  computed  as: 


tanA  ■  jiX'  cos 
d$ 


| 
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Substituting  the  expression  for  d\'/d$  from  (14)  gives  the  equation 
expressing  groundtrack  azimuth  on  the  sphere  as  a  function  of  the  orbit 
parameters. 


tanA  =  cosi  -  (p/P)cos2(t>  (16) 

(cos2<j>  -  cos2i)^ 

3.2  The  Mapping  Equations  For  The  Cylindrical  Projection 

Formulae  for  length  distortion  on  the  projection  surface  along  the 
meridians  (Mm)  and  along  the  parallels  (Mp)  are  defined  as: 


Mm  ■  _1  d^r 

(17a) 

R  d<|> 

Mo  *  1  dx 

(17b) 

R  cos$  dX 

3.2.1  The  mapping  equation  for  x 

To  derive  the  mapping  equation  for  x,  define  the  scale  to  be  true  in 
the  x  direction  at  4>j,  the  conformal  latitude. 

Mp  -  1 

dx  =  R  cos$ j  dX 

Integrating  with  the  constraint  that  x*0  when  X«0  gives  the  mapping 
equation  for  x. 

x  *  RXcos^j  (18) 

3.2.2  The  mapping  equation  for  v 

To  derive  the  parallel  spacing,  impose  the  condition  of  conformality 


» 


» 


f 


I 


< 


I 


I  < 


»  < 


I 


at  ip  j . 


Mm  *  Mp 

dy  ■  dx  (IS 

cos<frj/ dl\ 

W$/  $1 

Integrate  equation  (19)  and  substitute  from  (18)  replacing  X  with  X' 
to  make  y  a  linear  function  of  the  satellite  apparent  longitude. 
dX'/d<fr  is  a  constant  since  it  is  evaluated  at  <j>j. 


y  ”  RX "cos 

COS^j/ 


/dX'\ 
\d$  ft  i 


Because  conformality  is  defined  at  $j,  the  groundtrack  azimuth  at 
on  the  map  (dx/dy)  should  equal  the  azimuth  on  the  globe  (A). 
Rearranging  equation  (19),  again  using  X'  in  place  of  X  yields 


dx  =*/d_X  ^ 
dy  \d<j>  ft 


This  agrees  with  the  value  of  tanA  in  (15),  therefore 


(4T)  • 

\  d  4>  /4>  i 


COS$i 


Substituting  into  equation  (20)  gives  the  mapping  equation  for  y. 


y  *  RX'cosbi 
tanA, 


3.3  The  Mapping  Equations  For  The  Conic  Projection 


3.3.1  The  mapping  eauation  for  9 


Figure  2  shows  the  elements  of  the  projection  to  the  plane  in  the  (p,0) 


t 
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coordinate  system.  As  with  conventional  conics  in  the  normal  aspect, 
the  mapping  equation  for  0  is 

0  =  nX  (22a) 

In  the  derivation  of  p,  the  satellite  apparent  longitude  is  substituted 
for  earth-fixed  longitude. 

0  =  nX '  (22b) 

where 

n  is  the  constant  of  the  cone,  and 

0  is  the  angular  separation  of  the  meridians  on  the  projection 
surface. 

3-3.2  The  mapping  equation  for  p 

From  the  geometry  of  the  plane  triangle  in  Figure  2 
angle  C  -  180°  -  (0  +  S) 

where  S  is  the  azimuth  of  the  groundtrack  at  the  equator  on  the  pro¬ 
jection  surface.  Applying  the  sine  law  gives  a  relationship  between  p 
and  pQ,  the  radius  from  the  apex  of  the  cone  to  the  equator. 

P  *  PnSinS  (23) 

sin(0  4-  S) 

To  obtain  the  mapping  equation  for  p,  evaluate  the  length  pq  at  a 
conformal  latitude  and  substitute  into  equation  (23) . 

Formulae  for  length  distortion  in  the  (p,0)  system  are 


«»> 


Mm  ■  -_1  d£ 
R  d 


(24a) 
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(24b) 


The  minus  sign  in  Mm  indicates  that  p  decreases  as  $  increases. 


Differentiating  equation  (22b)  with  respect  to  p  and  X'  yields 


d9  = 

n 

(25) 

dX' 

d9  = 

n  dX' 

• 

d<p 

d 

Substituting  the 

value  of  dXVd<j>  from  equation  (15) 

d9  = 

n  tanA 

(26) 

• 

d<j> 

cos$ 

* 

Imposing  conformality  at  produces  an  expression  for  the  length  p^. 


Mm  =  Mp 

p  =  -cos<t>j/dX.A  /dg\ 

\de  )[<!*)*! 

Substituting  p  from  equation  (23)  and  dX'/d9  from  (25)  yields  ^ 

p  =  -cos^j  _1  (  Pn  sinS  \ 

n  d<t  \sin(9  +  S)/^ 

p.  =  cosjn  On  sinS  cos(91  +  S)  d  (9  +  S)  . 
n  sin‘-(91  +  S)  d(J> 

Substituting  d9/d$  from  equation  (26) 

» 

p  j  =  _ pn  sinS  tanAi _ 

sin(9j  +  S)  tan(9j  +  S) 

Comparing  this  with  equation  (23) 

I 

I 


Mp  =  p  d_9 
Rcosiji  dX 


tanAj  =  tan(9j  +  S) 
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Aj  =  0X  +  S  (27) 

Substituting  equation  (27)  into  (23)  evaluated  at  <j>j 

Pi  =  On  sinS  (28) 

sinAi 

Imposing  true  scale  at  <j>  produces  a  second  expression  for  pj. 

Mp  =  1 

p  =  Rcos^j  dA " 

de 

pj  =  Rcos^i  (29) 

n 

Substitute  equation  (29)  into  (28)  to  eliminate  pj  and  obtain  the 
equation  for  pQ. 


pQ  *  Rcosifri  sinAi 
n  sinS 


(30) 


Substitute  equation  (30)  into  (23)  to  eliminate  pQ  and  get  the 
mapping  equation  for  p. 

p  =  Rcosdn  sinAi  (31) 

n  sin (9  +  S) 

3.3.3  The  constant  of  the  cone  for  a  projection  with  two  confor¬ 
mal  latitudes 

Equation  (27)  gives  the  value  of  9  at  the  conformal  latitudes  ^  and 
$2  * 


91  -  A1  - 


S 
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®2  =  a2  ”  S 


Substituting  into  equation  (22b) 


9l  =  nXJ  =  Aj  -  S 


9o  —  nXo  Ao  S 


Subtracting: 


nX2  -  nX{  =  (A2  -  S)  -  (Aj  -  S) 

n  ■  a2  ,-jy.i 

Xf  -  Xi 


3. 3- 4  The  constant  of  the  cone  for  a  projection  with  one  conformal 
latitude 

As  $2  approaches  <j>j,  equation  (32)  becomes 


n  =  dA  =  (  cLa\  / dd>  \ 
dX'  \d4>j\dX-y 


Differentiating  equation  (16)  gives  the  expression  for  dA/d<f>. 


sec2A  dA  =  {(cos2<}>  -  cos2!)'5  (2 (p/P)sin<ji  cosij>) 
d<*> 

+  (cosi  -  (p/P)  (cos2(ji)  (4(cos2<()  -  cos2 i)-i4(-2sin<()  cos<ji)  } 

*  (C0S2<j>  -  COS2l) 


dA  =  sinift  cos#  { (p/P)  (cos2<fr  -  2cos2\)  4-  cost) 
dd>  (cos2<j>  -  cos2 1 )  3/2  (i  +  tan2A) 


Substituting  this  and  the  value  of  d<f>/dX'  from  (15)  into  the  equation 
for  n  evaluated  at  yields 


n  =  sinflj  cos4>]  {(p/P)  cos24>i  -  2cos20  +  cost)  cos<f>i 
(cos2<jii  -  cos2i)3^  (1  +  tanzAj)  tanAi 
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n  =  sin<i>T  cos24m  {(p/P)  cos24>i  -  2cos2i)  +  cosi} 

(cos2^  -  cos2t)3//2  (tanAj  +  tan3Aj) 

Expanding  the  denominator  using  the  expression  for  tanA  in  (16) 


(cos2(j>1  -  cos2i)3/2  /  cosi  -  (p/P)  cos2<fri  +  (cost  -  (p/P)  cos2^)3^ 

\  (COS24>2  -  COS2l)^  (C0S2(fl1  -  COS2l)3^2  J 


Thus 

n  =  sinjh  cos24n  {(p/P)  (cos24m  -  2cos2i)  +  cosi) _  (: 

(cos2^  -  cos‘!i)(cosi  -  (p/P)cos2tj>i)+(cosi  -  (p/P)cos2i j)  3 

The  equation  for  n  appearing  in  Snyder  (1981)  contains  additional 
expansions  of  the  denominator.  See  Appendix  A  for  a  continuation 
of  the  derivation. 


3.3.5  The  constant  of  the  cone  when  the  conformal  latitude  is  tht 
tracking  limit 

The  tracking  limit,  the  highest  latitude  reached  by  the  satellite 
equal  to  the  inclination  of  the  orbit.  To  determine  the  value  of 
replace  with  i  in  equation  (33). 

n  =  sini  cos2i  (cosi  -  (p/P)  cos2i) 

*  (cost  -  (p/P)  cos2i)3 

n  =  _ sin  i _  ( 

(1  -  (p/P)  cosi)2 

3.3.6  The  circle  of  groundtrack  tangencv 

If  the  projection  is  to  be  constructed  manually,  the  plotting  of 
groundtracks  is  made  easier  by  computing  the  circle  of  tangency. 
is  a  circle  of  radius  pg  to  which  all  groundtracks  are  tangent. 
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1 

* 
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From  the  geometry  of  the  plane  triangle  in  Figure  2 

> 

ps  =  pQ  sinS 

Substituting  the  value  of  pQ  from  equation  (30) 

pg  =  Rcosjn  sinAi  (35) 

n 

The  direction  of  the  groundtracks  is  obtained  from  the  sign  of  either 
Ps  in  equation  (35)  or  S  in  equation  (27).  A  negative  sign  indicates 
a  west  azimuth  when  the  satellite  is  travelling  northward,  or  east 

azimuth  when  it's  going  southward.  ' 

Equation  (35)  is  needed  to  compute  pg  because  equation  (31)  breaks 
down  when  <j>  equals  the  tracking  limit.  Because  parallel  spacing  was 

I 

derived  based  on  the  satellite  groundtrack,  radii  poleward  of  the 
tracking  limit  are  undefined,  i.e.,  the  sine  of  a  computed  in 
equation  (2)  is  greater  than  unity. 


4.  DISTORTION  ANALYSIS 


The  Satellite  Tracking  projections  were  developed  by  defining  the 
relationship  between  parameters  on  the  projection  surface  and  para¬ 
meters  on  the  sphere  at  specific  locations,  i.e.,  conformality  and 
true  scale  at  one  or  two  latitudes.  Now  that  the  mapping  equations 
have  been  developed,  general  equations  describing  the  distortion 
properties  of  the  projections  can  be  computed  from  fundamental  quan¬ 
tities  developed  using  the  following  transformation  matrix. 

E 

F 

G 

where 

u,v  are  the  parametric  quantities  on  the  projection  surface, 
either  (x,y)  or  (p,8), 

<J>,A  are  the  parameters  on  the  sphere, 

E',  F",  G'  are  the  fundamental  quantities  of  u,v  with 
respect  to  the  projection  surface, 

E,  F,  G  are  the  fundamental  quantities  of  u,v  with  respect 
to  the  sphere. 

Also  used  in  the  distortion  analysis  are  e,  f,  and  g,  the  fundamental 
quantities  of  $,A  with  respect  to  the  sphere. 


du2 

2  du  dv 

dv2 

d<|>2 

d<j>  dif> 

d<j>2 

du  du 

du 

dv  +  du  dv 

dv  dv 

d$  dX 

d4> 

dX  dX  d(J> 

dX  d<f> 

du2 

2  du  dv 

dv2 

dX2 

dX  dX 

dX2 
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The  following  fundamental  quantities  for  the  sphere  and  plane  (pro¬ 
jection  surface)  are  known  through  the  equations  defining  differen¬ 
tial  length 


4>,X  on  the  sphere:  ds2  =  R2 (d^*) 2  +  R2cos2(ji  (dx)2 

e  =  R2  f  —  0  g  =  R2cos2<ji 

x,y  on  the  plane:  dS2  =  (dx)2  +  (dy)2 

E'  =  1  F'  =  0  G'  -  1 


p,0  on  the  plane:  dS2  *  (dp)2  +  p2(d0)2 

E'  -  1  F'  =*  0  G'  =  p2 


4.1  Cylindrical  Projection  Distortion  Equations 

4.1.1  The  first  fundamental  quantities 

Recalling  equations  (15),  (18),  and  (21),  and  evaluating  their  partials 


dX'  =  tanA  (15) 

d$  cos((> 

x  =  RXcos4>j  (18) 

y  =  RX^cosifri  (21) 

tanAj 

dx  =  0 
dtp 

dx  =  Rcosif! 
dX 

dy  ■  Rcosdn  dX'  =  Rcosiln  tanA 
d<J>  tanAj  dtp  tanAj  cosij> 


dy  =  0 
dX 
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Substitute  the  known  values  into  the  transformation  matrix. 


E 

0 

0 

R2cos24n  tan2A 

1 

tan2Aj  cos2$ 

F 

_ 

0 

R2cos2d> 

i  tanA 

0 

0 

tanA  i 

cos<t> 

G 

R2cos24>i 

0 

0 

1 

- 

-  - 

E  = 

R2cos24ii  tan2A 
tan2Aj  cos^4> 

(36a) 

F  = 

0 

(36b) 

G  * 

R2COS2<jl  J 

(36c) 

These  are  the  fundamental  quantities  of  the  cylindrical  projection. 


4.1.2  Scale  distortion  in  the  meridians  and  parallels 


Mm 


R2cos2Jul 


tan2Ai 


cos 


Mm  *  cosifri  tanA 
tanA1  cos$ 


Mp  = 


(37) 


Mp  =  cos<t>i 
cos$ 


(38) 


4.1.3  Area  distortion 


Da  -  /EG  -  ?2\h 
\eg  -  f2/ 


{- 


Rucos4d>i  tan2A\fe 
tan2Ai  cos  2 4 

- RW* 


Da  =  cos2<fri  tanA 
tanAi  cos2$ 

4.1.4  Angular  distortion 

A  derivation  of  the  u>  to  Q  equation  is  given  in  Appendix  B. 


cos  ft  *  i/?  cos  <»  cotta 

(E  cos2f  cot2u  +  0*5 


Rcostfri  tanA  cos»  cotta 

*  _  tanAi  costfr _ 

/R2cos24n  tanzA  cos^tfe  cot^oj  +  R2cos2(frA  H 
l  tan2Aj  cos2$  J 

cosft  *  tanA  cotta 

(tan2A  cotzu>  +  tan^Aj)^ 


4.2  Conic  Projection  Distortion  Equations 
4.2.1  The  first  fundamental  quantities 

Recalling  equations  (15),  (22a),  (22b),  and  (31)  and  evaluating 
their  partials: 

dl*  ■  tanA 
d$  cos<$ 

9  =  nX 

9  *  nA " 
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0  =  Rcosifti  sinAi  (31) 

n  sin(9  +  S) 

de  =  o 

d<j> 

d6  =  n 
d  A 

dp  =  Rcosdn  sinAi  d  f _ 1 _ \ 

d  4>  n  d<j>ysin(  9  +  S  )J 

=  -Rcosifri  sinAi  cos (9  4-  S)  d  (9  +  S) 
n  sinz( 9  +  S)  d<J> 

Because  the  mapping  equation  for  p  was  developed  using  satellite 
apparent  longitude,  the  value  for  0  in  dp/ d d>  is  gotten  from  equation 
(22b).  Substituting  from  (22b)  and  (15)  yields 


d  (9  +  S)  =  d_(n A 0  +  d£5  *  n  d_A'  =  n  tanA 
d  dij)  d  4>  d  d>  cos<j> 


dp  =  -Rcos<fri  sinAi  tanA 
d d>  cos<j>  sin(9  +  S)  tan(8  +  S) 

dp  =  0 
dA 


Substituting  the  partials  into  the  transformation  matrix 


E 

R2cos2<)>t  sin2A1  tan2A  0 

0 

1 

cosz<J>  sinz(0  +  S)  tanz(9  +  S) 

F 

* 

-n  Rcosifei  sinAi  tanA 

0  cosip  sin(9  +  S)  tan(0  +  S) 

0 

0 

G 

0  0 

n2 

p2 

—  m 

E  = 


_ R2 

cos24> 


cos2^  sin2A]  tan 2 A 
sin*(9  +  S)  tan2(6  + 


S) 


(41a) 
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F  -  0  (41b) 

G  =  n2p2  (41c) 

Since  p  is  known  as  a  function  of  substitute  equation  (31)  into 
(41c)  to  obtain  the  equation  for  G. 


G  =  R2cos2ifri  sin2Ai 
sin2 (9  +  S) 

These  are  the  fundamental  quantities  of  the  conic  projection. 


(4  Id) 


4.2.2  Scale  distortion  in  the  meridians  and  parallels 
Mm  =  /E\J5 


/  R2cos2Ai  sin2Ai  tan2A _ \  % 

cos 2d>  sln*(9  +  S)  tan2(9  +  S) 

/ 


Mm  -  _ cosifti  sinAi  tanA 

cos<Ji  sin(9  +  S)  tan(9  +  S) 


(42) 


Mp  =  (G\h 


(f 


(r^os2^  sin2Ai^  h 
s^-n2(8  +  S) 


R2cos2<t> 


Mp  =  cosj^  sinAi 
cos<j>  sin(9  +  S) 


(43) 


4.2.3  Area  distortion 


Da 


(eg-fAh 

\eg  -  f2/ 


Da  = 


4.2.4 

cosft 


cosfl 


R  cosH<fri  sin^Ai  tan2A 
cosz(j>  sin^CQ  +  S)  tan2Q  +  ¥) 
R^cos2* 


cosit, 

5  ^  r+i  o  ■(  n 


sin2Ai  tanA 


cos2(ji  sinz(9  +  S)  tan(9  +  S) 


Angular  distortion 

_ /E  COS<j>  cotco 

(E  cos2<j>  cot2co  +  G)^ 


£ 


Rcos<j>i  sinAT  tanA  cos4>  cotco 
costfe  sin (9  +  tan (9  +  S) 


sin-^Ai 

cosz<j>  sin2  (9  +  S)  tan2 (9  +  S) 


sin^Ai 


cos“ 
Zi 


LL 


tan2A 


cos 

rz 


-± 


cot2!]) 


+  R2 


COS2 1^  T 


sin2 (0  +  S) 


T 


tanA  coto)  _  , 

CtanzA  cot2(j  +  tan2  (9  +  S))  * 


5.  RESULTS  AND  ANALYSIS 


Inspecting  the  distortion  equations  for  both  projections  reveals 
that  they  do  satisfy  the  conditions  of  true  scale  and  conformality 
imposed  during  the  derivation  of  the  mapping  equations. 

<P1 

Mm*  =  Mp . 

*2 

Mm.  =  1 
<J>1 

* 1 

ft .  =  w  . 

♦  l  $  i 

ft.  =  W  4. 

*2  $2 

As  a  further  check.  Da  =*  Mm  Mp. 

The  following  analysis  is  illustrated  with  plots  of  the  various  distor¬ 
tions  for  the  Landsat  orbit.  Tables  containing  the  numerical  distor¬ 
tion  values  appear  in  Appendix  C. 

The  parameters  used  in  the  calculations  are 

satellite  revolution  period  =  103.267  minutes, 
earth  rotation  period  =  1440.0  minutes, 
satellite  inclination  =  99?092. 

5. 1  The  Cylindrical  Projection 

The  meridian  scale  distortion,  Figure  3,  is  a  function  of  the  secant  of 
the  latitude.  However,  because  the  spacing  of  the  parallels  is  based  on 
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Figure  3.  Meridian  length  distortion  at  standard  parallels  of  0°,  40°,  60°,  and  80 


the  satellite  orbit,  the  effect  of  the  orbit  appears  in  the  tanA  term. 
Since  the  azimuth  of  the  groundtrack  on  the  sphere  increases  with 
latitude,  the  large  range  of  Mra  reflects  the  increased  parallel  spacing 
required  to  introduce  sufficient  angular  distortion  to  keep  the  projected 
azimuth  straight. 

The  <f>j  subscripted  terms  in  all  of  the  distortion  equations  allow  the 
scale  factor  to  equal  unity  at  the  standard  parallel. 

The  parallel  scale  distortion,  Figure  4,  is  a  function  only  of  the  secant 
of  the  latitude.  This  indicates  that  latitude  length  is  independent 
of  its  projected  distance  from  the  equator  because  the  meridians  are 
cast  as  equidistant  vertical  lines.  Area  distortion,  the  product  of 
Mm  and  Mp,  is  shown  in  Figure  5. 

Figures  6,  7,  and  8  show  the  effect  of  angular  distortion  at  standard 
parallels  of  20°,  45°,  and  70°.  Conformality  insures  that  the  azimuth 
on  the  map  equals  the  azimuth  on  the  sphere  at  the  standard  parallel. 
Because  azimuth  is  a  constant  on  the  projection,  but  increases  with 
latitude  on  the  sphere,  constraining  the  groundtracks  as  straight  lines 
causes  map  azimuth  to  be  greater  than  sphere  azimuth  below  the  standard 
parallel  and  less  than  sphere  azimuth  above  it.  This  effect  can  also 
be  seen  in  the  scale  distortions:  Mm  is  less  than  Mp  below  ^  and 
greater  than  Mp  above  <j>j  due  to  the  relative  sizes  of  A  and  Aj  in 
equation  (37). 


at  standard  parallels  of  0°,  40  y  60  t  and  80 
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5.2  The  Conic  Projection 

The  equations  for  scale  distortion  are  similar  to  those  for  the 
cylindrical  projection;  but  here,  the  azimuth  of  the  groundtrack 
changes  on  the  projection  as  well  as  on  the  sphere.  This  is  reflected 
in  the  addition  of  the  sinAj/sin(9  +  S)  term  in  both  equations,  which 
relates  the  azimuth  at  the  standard  parallel,  Aj,  to  its  general  value 
of  (9  +  S).  Plots  of  the  meridian  and  parallel  scale  distortions  for 
the  tangent  conic  projection  are  shown  in  Figures  9  and  10;  area  dis¬ 
tortion  is  shown  in  Figure  11. 

Inspecting  the  tables  of  the  scale  distortion  values  in  Appendix  C 
reveals  that,  for  the  tangent  cone,  Ifai  is  greater  than  Mp  everywhere 
except  at  the  standard  parallel  where  they  are  equal.  For  this  to  be 
true,  comparing  equations  (42)  and  (43),  the  azimuth  on  the  sphere.  A, 
must  always  be  greater  than  or  equal  to  the  azimuth  on  the  plane, 

(9  +  S).  This  is  confirmed  in  the  angular  distortion  plots  in  Figures 


12,  13,  and  14. 
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Figure  13.  Conic  projection  azimuth  distortion. 

Standard  parallel  =  45  degrees. 
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6.  CONCLUSIONS 


The  purpose  of  this  thesis  is  to  derive  general  distortion  formulae 
for  the  Satellite  Tracking  map  projections.  The  following  conclu¬ 
sions  can  be  drawn: 

1.  Equations  for  length,  area,  and  azimuth  distortion  for  the 
Satellite  Tracking  projections  were  derived. 

2.  Limitations  in  the  use  of  these  projections  are  now  quanti¬ 
tatively  defined,  e.g.,  the  size,  shape,  and  orientation  of  a 
Landsat  scene. 

3.  The  cylindrical  projection  is  preferred  when  the  entire  earth 
is  to  be  shown  because  it  allows  true  scale  to  be  defined  both 
above  and  below  the  equator. 

4.  The  conic  projection  is  preferred  for  depicting  small  areas 
because  conformality  can  be  defined  at  any  two  parallels 
within  the  area  of  interest. 

Recommendations  for  further  study: 

These  projections  are  valid  only  under  the  assumption  of  a  spherical 
earth  and  circular  satellite  orbit. 

1.  What  errors  in  the  plotted  groundtrack  are  introduced  as  a 
result  of  these  assumptions. 

2.  Can  the  projections  be  modified  to  accept  a  non-circular  orbit. 
Can  the  transformation  be  made  from  an  ellipsoidal  earth. 


3. 
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Appendix  A.  A  continuation  of  the  derivation  of  the  constant  of  the 
cone  for  a  projection  with  one  conformal  latitude. 

Recalling  equation  (33): 

n  =  sindn  cos2$i  { (p/P) (cos2dn  -  2cos^i)  +  cost) _ 

(cos^j  -  cos^ i )  (cos t  -  (p/P)cos2<j>j)  +  (cos i  -  (p/Pjcos2^)  3 

Expanding  the  denominator: 

(cos2^  -  cos2i)(cosi  -  (p/P)cos2<j>i)  + 

(cos3!  -  3(p/P)coS2({)1  cos2!  +  3(p/p)2  COS4d>l  COS!  -  (p/P)3  COS6<p:) 

=  -2  (p/P)  COS2<()  1  cos2!  +  3(p/P)2COS4'j>1  COS!  - 
(p/p)3  cos^ijlj  4-  COS2<))2  COS!  -  (p/P)coS4(j>i 

=  COS2(()1{-2  (p/P)cOS2!  +  3(p/p)2  COS2{J)  |  COS!  - 

(p/P)  3  COS4^  +  COS!  -  (p/p)  COS2(j)  |  } 

=  COS2({(1  (~(p/P)cOS2<jl1  {(p/P)2  COS2(J>1  -  2  (p/P )  COS !  +  1}  + 

COS!  {(p/P)2  COS2^!  -  2 (p/P) COS !  +  1}) 

=  COS2^  ((COS!  -  (p/P)cOS2<J>j){(p/P)2  cos2<j)1  -  2  (p/P)  COS  1  +  1}) 

=  COS2  j  ((COS!  -  (p/P)cOS2<Jl1){(p/P)  ( (p/p)  COS2  d>  J  -  2coS!)  +  1}) 

Substituting  this  back  into  equation  (33)  yields  the  final  form  of  n 

n  =  _ sin<ti  {(p/P)  (cos2<tn  -  2cos2!)  +  cos!) 

(cosi  -  (p/P)cos2(jii)  { (p/P)  ( (p/P)cos2  jij  -  2cosi)  +  1} 

This  form  is  useful  to  test  the  condition  where  =  90°  and  the 
projection  becomes  azimuthal.  This  could  not  be  done  using  equation 
(33)  because  of  the  cos2<t>.  terra  in  the  numerator. 
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Appendix  B.  The  Derivation  of  the  Angular  Distortion  Equation 


On  tie  projection  surface: 

cosft  =  v*E  _d<£ 
dS 

=  v'E  dd> 

A  (<4H  +  G  (dA)z 

cosfi  =  /E~  d 

/E7d*p+  G  dX 
\dx) 


On  the  sphere: 


ds  cosw  =  R  d$> 

ds  sinw  =  R  cos<£  dX 

d^>  =  R  cos<t)  costii 
dA  R  sinw 

d<fr  =  cos<t>  cotw 
dA 

Substituting: 

cosft  =  /iE  cos<t>  cotu 

/E  costii  cotzuj  +  G 
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Appendix  C.  Distortion  Tables. 


I 

Table  1 

I 

Table  2 


Table  3 


Cylindrical  projection  length  and  area  distortion  for  th 
Landsat  orbit. 


Conic  projection  length  and  area  distortion  for  the 
Landsat  orbit. 


Cylindrical  projection  azimuth  distortion  for  the  Landsa 
orbit. 


Table  4.  Conic  projection  azimuth  distortion  for  the  Landsat  orbi 
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Table  1.  Cylindrical  projection  length  and  area  distortion  for  the 
Landsat  orbit. 

•  < 


Standard 

Parallel 

0° 

o 

O 

60° 

00 

O 

o 

Meridian 

Length  Distortion 

0° 

1.0 

.66762 

.31363 

.01816 

10 

1.02179 

.68217 

.32047 

.01855 

20 

1.09298 

. 72969 

.34279 

.01985 

30 

1.23456 

.82421 

.  38720 

.02242 

40 

1.49787 

1.0 

.46978 

.02720 

50 

2.01389 

1.34451 

.63162 

.03657 

60 

3.  18846 

2. 12866 

1.0 

.05790 

70 

6.89443 

4.60283 

2. 16231 

. 12519 

80 

55.07144 

36.76658 

17.27213 

1.0 

80.908 

00 

ao 

ao 

00 

Parallel  Length  Distortion 


0 

1.0 

. 76604 

.5 

. 17365 

10 

1.01543 

.77786 

.50771 

.17633 

20 

1.06418 

.81521 

.53209 

. 18479 

30 

1.15470 

.88455 

.57735 

.20051 

40 

1.30541 

1.0 

.65270 

.22668 

50 

1.55572 

1. 19175 

.77786 

.27015 

60 

2.0 

1.53209 

1.0 

.34730 

70 

2.92380 

2.23976 

1.46190 

.50771 

80 

5. 75877 

4.41147 

2.87939 

1.0 

80.908 

6.32830 

4.84776 

3. 16415 

1.09890 

0 

1.0 

10 

1.03756 

20 

1. 16212 

30 

1.42555 

40 

1.95533 

50 

3. 13306 

60 

6.37691 

70 

20. 15798 

80 

317. 14380 

80.908 

on 

Distortion 


.51142 

. 15682 

.53063 

. 16270 

.59485 

. 18240 

. 72906 

.22355 

1.0 

.30663 

1.60232 

.49131 

3.26130 

1.0 

10.30926 

3.  16109 

162. 19480 

49. 73312 

on 


»  < 


.00315 

.00327 

.00367 

.00449 

.00617 

.00988 

.02011 

.06356 

1.0 


» 


00 


50 


»  4 


Table  2.  Conic  projection  length  and  area  distortion  for  the 
Landsat  orbit. 


I  < 


Standard 

Parallel 

10° 

30° 

50° 

O 

O 

» 

Meridian  Length  Distortion 

0° 

1.00313 

1.07991 

1. 71646 

23.26322 

10 

1.0 

1.00398 

1.35754 

6.92079 

» 

20 

1.04275 

.97729 

1. 14552 

3.37126 

30 

1.  14572 

1.0 

1.02572 

2.06490 

40 

1.34649 

1.08700 

.97692 

1.46129 

50 

1.73949 

1.28047 

1.0 

1. 15710 

60 

2.60194 

1.69556 

1.12369 

1.01679 

70 

5.08699 

2.73190 

1.44871 

1.0 

> 

80 

29.00155 

9.41190 

3.09745 

1.22864 

80.908 

03 

<30 

00 

00 

( 


< 


I 


Parallel  Length  Distortion 

0 

.99671 

.99621 

1. 17680 

4.05740 

10 

1.0 

.96599 

1.05324 

2.22532 

20 

1.03510 

.96693 

.98251 

1.57647 

30 

1.10816 

1.0 

.95174 

1.26294 

40 

1.23355 

1.07194 

.95670 

1.09489 

50 

1.44188 

1.19859 

1.0 

1.00824 

60 

1.80333 

1.41446 

1.09269 

.97861 

70 

2.51103 

1.79895 

1.25919 

1.0 

80 

4.20804 

2.42253 

1.47016 

1.06676 

80.908 

4.23325 

2.24930 

1.35775 

1.04704 

Area  Distortion 

0 

.99983 

1.07582 

2.01993 

94.38810 

10 

1.0 

.96984 

1.42982 

15.40095 

20 

1.07935 

.94497 

1. 12548 

5.31469 

30 

1.26964 

1.0 

.97622 

2.60784 

40 

1.66097 

1. 16520 

.93462 

1.59995 

50 

2.50814 

1.53477 

1.0 

1.16663 

60 

4.69216 

2.39830 

1.22784 

.99504 

70 

12.77360 

4.91455 

1.82420 

1.0 

80 

122.03977 

22.80060 

4.55374 

1.31066 

80.908 

00 

00 

00 

oo 

i 
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Table  3-  Cylindrical  projection  azimuth  distortion  for  the  Landsat 
orbit. 


Latitude 

0° 

O 

O 

60° 

O 

O 

oo 

Sphere 

Azimuth 

o 

Map  Azimuth 

Standard  Parallel 

0 

=  20° 

O 

0° 

0.0 

0.0 

0.0 

0.0 

10 

10.26495 

8. 96897 

6.48077 

1.08436 

20 

20.49676 

18.04498 

13. 19640 

2.23744 

30 

30.66688 

27.32916 

20.40258 

3.54643 

40 

40.75502 

36.90925 

28.39457 

5. 14695 

50 

50.75153 

46.84937 

37.51585 

7.29022 

60 

60.65803 

57. 17734 

48.13385 

10.53262 

70 

70.48670 

67.87198 

60.53498 

16.43235 

80 

80.25838 

78.85585 

74.69298 

31.33252 

90 

90.0 

90.0 

90.0 

90.0 

Standard  Parallel 

II 

Ln 

O 

0 

0.0 

0.0 

0.0 

0.0 

10 

12.03481 

10.52545 

7.61675 

1.27646 

20 

23.75263 

20.98287 

15.43133 

2.63343 

30 

34.91723 

31.31478 

23.64675 

4. 17284 

40 

45.41324 

41.48222 

32.47151 

6.05278 

50 

55.23919 

51.46883 

42.10802 

8.56439 

60 

64.47477 

61.28085 

52. 71885 

12.34601 

70 

73.24646 

70.94427 

64.36273 

19. 14672 

80 

81. 70266 

80.50024 

76.91131 

35.62808 

90 

90.0 

90.0 

90.0 

90.0 

Standard  Parallel 

o 

O 

0 

0.0 

0.0 

0.0 

0.0 

10 

22.57733 

19.91899 

14.61773 

2.48838 

20 

40.63871 

36.79633 

28.29625 

5. 12597 

30 

53.70224 

49.87562 

40.49658 

8.09850 

40 

63. 18850 

59.89055 

51. 14122 

11.68437 

50 

70.41251 

67. 78976 

60.43412 

16.36862 

60 

76.24248 

74.30808 

68.67779 

23. 11683 

70 

81.22568 

79.95673 

76. 17587 

34. 10376 

80 

85.72366 

85.09605 

83.20187 

54.41874 

90 

90.0 

90.0 

90.0 

90.0 

52 


Table  4.  Conic  projection  azimuth  distortion  for  the  Landsat  orbit. 


Latitude _ 0^ _ 40^ _ 60^ _ 80° 

Map  Azimuth 

Sphere 

Azimuth  Standard  Parallel  =  20° 


0° 

0?0 

0?0 

0?0 

o!o 

10 

9. 71601 

9.55806 

7.64103 

1.95248 

20 

19.46493 

19. 16619 

15.47873 

4.02519 

30 

29.27632 

28.86987 

23.71462 

6.36914 

40 

39. 17314 

38.70529 

32.55515 

9.21472 

50 

49. 16891 

48.69490 

42. 19983 

12.97502 

60 

59.26561 

58.84354 

52.80775 

18.51405 

70 

69.45275 

69.13623 

64.43463 

27.97665 

80 

79.70805 

79.53849 

76.95195 

47.63448 

90 

90.0 

90.0 

90.0 

90.0 

Standard  Parallel  * 

45° 

0 

0.0 

0.0 

0.0 

0.0 

10 

7. 70744 

9.95383 

9.44635 

4.08611 

20 

15.60828 

19.91319 

18.95443 

8.38833 

30 

23.89994 

29.88297 

28.58072 

13. 16523 

40 

32. 78316 

39.86682 

38.37095 

18.77562 

50 

42.44963 

49.86671 

48.35445 

25.77252 

60 

53.04910 

59.88269 

58.53894 

35.05824 

70 

64.62944 

69.91288 

68.90689 

48.06408 

80 

77.06188 

79.95362 

79.41530 

66.48014 

90 

90.0 

90.0 

90.0 

90.0 

Standard  Parallel  = 

'vl 

o 

o 

0 

0.0 

0.0 

0.0 

0.0 

10 

1.76150 

7.52606 

9.63174 

8.70415 

20 

3.63232 

15.25417 

19.30566 

17.53739 

30 

5.75013 

23.39270 

29.05984 

26.62384 

40 

8.32610 

32.15784 

38.92427 

36.07508 

50 

11.74211 

41.76285 

48.91710 

45.97808 

60 

16.80910 

52.38369 

59.04168 

56.37768 

70 

25.60347 

64.09087 

69.28501 

67.25653 

80 

44.68725 

76.75732 

79.61826 

78.52028 

90 

90.0 

90.0 

90.0 

90.0 

Appendix  D.  FORTRAN  program  and  sample  output. 

This  program  computes  values  for  constructing  the  graticules  of  the 
Satellite  Tracking  projections.  The  type  of  the  projection  is  deter¬ 
mined  by  the  values  in  the  assignment  statements  in  the  "initialize 
variables"  section.  For  example,  if  PHI1  equals  zero,  or  PHI1  equals 
-PHI2,  a  cylindrical  projection  is  produced;  otherwise  the  result  is  a 
conic. 

Output  is  included  for  cylindrical  and  tangent  conic  projections  with  a 
standard  parallel  of  30°.  Plots  of  the  graticules  are  shown  in  Figures 
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CYLINDRICRL  PROJECTION 
STANDARD  PARALLEL  =  30° 


Figure  15.  Cylindrical  projection  graticule. 
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